Phase diagram of the t-U-J\-Ji chain at half filling 
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We investigate the half-filled Hubbard chain with additional nearest- and next-nearest-neighbor 
spin exchange, Ji and J2, using bosonization and the density-matrix renormalization group. For 
J2 = we find a spin-density-wave phase for all positive values of the Hubbard interaction U and 
the Heisenberg exchange J\. A frustrating spin exchange J2 induces a bond-order- wave phase. For 
some values of Ji, J2 and U, we observe a spin-gapped metallic Luther-Emery phase. 
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I. INTRODUCTION 



The Hubbard chain is the archetype of one-dimensional 
strongly correlated electron systems. At half band- 
filling and for all values of the Hubbard interaction U, 
it exhibits insulating spin-density-wave (SDW) behavior, 
marked by a critical behavior of the spin correlations. In 
a weak-coupling picture, this insulating behavior is gener- 
ated by umklapp scattering, while in strong coupling, the 
opening of the Mott-Hubbard gap leads to behavior of the 
spin degrees of freedom governed by an effective Heisen- 
berg chain. These perturbative results are reinforced by 
the exact Bethe-Ansatz solution.— Hubbard-type models 
are relevant to a wide variety of one-dimensional materi- 
als, such as polymers^ cuprates^ or TTF-TCNQ^ 

One important experimental question is to what ex- 
tent the spin correlations remain critical when addi- 
tional interactions are present. It is well known that 
any dimcrization or sufficiently large frustration can lead 
to a spin gap in the Heisenberg- type spin models. A 
nearest-neighbor Coulomb repulsion^^ an alternating lo- 
cal potential^ or a second-neighbor hopping^^ can 
lead to a spin gap in models for itinerant interacting elec- 
trons. 

The Hubbard model with a nearest-neighbor antiferro- 
magnetic exchange in two dimensions is of interest in the 
context of the high-T c cuprates. In particular, spin-liquid 
states^ and gossamer superconductivity^ at and near 
half-filling have been proposed as necessary precursors 
to high-temperature superconductivity at higher doping. 
Since it is not clear whether such states are present in 
sufficient strength and for sufficiently wide parameter 
regimes in the pure Hubbard or t-J models, additional 
interactions, including a spin exchange, have been pro- 
posed to be relevant^ 

In this work, we investigate the effect of two additional 
terms on the phase diagram of the half-filled (average 
electron occupation (n) = 1) Hubbard chain, namely 
explicit antifcrromagnctic exchange interactions between 
nearest neighbors and between next-nearest neighbors. 



The Hamiltonian is given by 

H = -tJ2 { c L c i,a + h - c ) + u^2 n ^i 

i,a i 

+ J\ SiSi+i + J2 SiSi+2, (1) 

i i 

where cj a (c; i0 -) creates (annihilates) an electron with 
spin a at site i, n.i^ = c\ a c ia , and Si is the spin op- 
erator on site i: S? = \ £ CT ^ 4^>V _ The indcx 
a = x,y,z, and a^ a , are the Pauli matrices. Here t is 
the hopping amplitude and U the strength of the on-site 
Coulomb interaction. The antifcrromagnctic Heisenberg 
parameters J\ and Ji correspond to nearest- and next- 
nearest-neighbor exchange, respectively. 

The unfrustrated ( J2 = 0) version of this model has 
previously been investigated both analytically and nu- 
merically. In particular, a generalized model with an 
anisotropic Heisenberg coupling was investigated in Ref. 

using bosonization. While this work concentrated pri- 
marily on the case of ferromagnetic exchange, isotropic 
antifcrromagnctic exchange was included in a phase that 
is marked as 'dimcr long-range order', which corresponds 
to a bond-order wave (BOW) in our notation; see below. 
The phase diagram from bosonization of the isotropic 
antiferromagnetic exchange was considered explicitly in 
Refs. [HI and fTsl. supported by numerical calculations us- 
ing the transfer-matrix renormalization group (TMRG)i^ 
and exact diagonalization^ The phase diagram found 
contains two phases: a bond charge-density-wave phase 
(our BOW phase) at sufficiently small U for all 3\ , and 
a SDW at larger U. The critical value of U c goes to 
zero at small and large J\ and reaches a maximum value 
U c /t sa 0.35 at intermediate J\. 

In this work, we reexamine the bosonization treatment 
of the t-U-Ji model in the weak-coupling regime, includ- 
ing the renormalization of the coupling constants within 
the mean-field approximation. In addition, we consider 
the effect of the additional frustrating exchange J2 , which 
allows us to explicitly induce the bond-order phase and 
to make contact with the known phase diagram of the 
frustrated Heisenberg chain at large U. We also carry 
out high-precision ground-state density-matrix rcnormal- 
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ization group (DMRG) calculations, which allows us to 
explore the phase diagram numerically exactly. Both the 
revised bosonization and the DMRG calculations indi- 
cate that a BOW phase is not present for J 2 = 0; the 
system is in a SDW phase for all positive J\ and U . We 
show that a BOW phase can be induced by turning on 
J 2 positively, with the critical value required depending 
on U and J\. At larger values of J 2 , we find additional 
phases, including a spin-gapped metallic phase which we 
identify as a Luther-Emery phase. 

The paper is organized as follows: In Sec. [Til we dis- 
cuss the bosonization calculation and the resulting phase 
diagram. Sec. IIIII contains our numerical DMRG results 
and compares and contrasts the behavior obtained with 
that predicted by bosonization. In Sec. IIVI we discuss 
the overall phase diagram of the model in terms of the 
results from the two methods as well as the implications 
of our findings. 



II. FIELD THEORY 

We start our investigation with an analytical treat- 
ment of our model for small couplings, U, Ji, J 2 "C t. For 
simplicity, we take K = 1 everywhere. 



A. Linearization of the spectrum 

In terms of fermion operators the Hamiltonian (UJ) has 
the form 



H.c 



1=1 i 



' c i,t c i,t c i+£,i c i+e,i 



c i,l c i,l c i+i,l c i+i,l 
' C ll C J,l C i+f,T C i+f,T 



(2) 



For low temperatures and for excitations at low energies, 
it is enough to consider a restricted Hilbert space which 
contains only states close to the Fermi surface. In one 
dimension, the Fermi surface consists only of two points, 
k = ±kp. Around the Fermi points, the spectrum can be 
linearized and one can introduce left-moving and right- 
moving fermions corresponding to the states near — fcp 
and +A:p, respectively, 



Ci+^.+eM^+AO + -iMft+AO (3 ) 



for I = 0, 1,2. Here E4 is the coordinate vector of the 
site i and a is the lattice constant. For the half- filled sys- 
tem, &f = 7r/2a. Therefore, the left- and right-moving 
fermions have the phase factor e ±l£7r / 2 , for different val- 
ues of t. When written in terms of the chiral fermions 
c i+e.a,±, each interaction term of Hamiltonian ^ splits 



into four new terms. Two of them correspond to forward- 
scattering processes whose couplings are denoted by fl 2 
and 514 in standard g-ology notation.— In addition, there 
are two backward-scattering processes which describe 
"true" backward scattering (<?i-processes) and umklapp 
scattering ((^-processes). Due to the SU(2) symmetry of 
the spin sector, all processes depend only on the relative 
spins of the scattering electrons. This is denoted by the 
subscripts || and _L if the scattering electrons have the 
same or opposite spins, respectively. The relation be- 
tween the fl-ology parameters and the couplings of our 
original model is 



.92 _L 
53_L 
941. 



U 
U 
U 
U 



Ji/2 - 3J 2 /2, 
Ji/2-3J a /2, 
3Ji/2 - 3J 2 /2, 
3Ji/2 - 3J 2 /2, 



(4a) 
(4b) 
(4c) 
(4d) 



and 



9l\\ = ~ Ji/2+ J 2 /2, (5a) 

Sal, =Ji/2 + J 2 /2, (5b) 

ff3 |l = - J x /2 + J 2 /2, (5c) 

fl 4 || =Ji/2 + J 2 /2. (5d) 

In order to analyze the low-energy fl-ology model, we 
apply the bosonization method. 

B. Bosonization of the Hamiltonian 

First, we introduce the continuous chiral fermion 
fields ipa-,±(x) by making the replacement Ci^ a ,±/Va — > 
ip<r.±{x). The bosonization of the on-site interaction is 
straightforward. Using Abclian bosonization, we intro- 
duce the chiral boson phase fields 4> a ,±(x) via 

Vv,±(z) = -i=F ± e ±i2 *-±<*), (6) 



where F± are the so-called Klein factors which ensure the 
anti-commutation relations of the fermion fields. The 
symmetric and antisymmetric combination of the spin- 
dependent boson fields, 4>c,± = 4>t,± + 4>i,± an d 4>s,± = 
<t>1 t ± — 0x,±, correspond to the collective charge and spin 
modes, respectively. 

In order to bosonize the non-local processes, one must 
expand the fermion fields with respect to the lattice con- 
stant. The bosonized form of the fl-ology Hamiltonian 
density, up to leading order in the expansion with re- 
spect to the lattice constant, is 

H(0) ( X )=^-J2 [M^crf +V cr (d x( f> s , r ) 2 ] 



2 71 

9p 



r=± 



+ i^{d x 4>c, + ){d x 4>c,-) - cos(2^) 
^ ^(d x <p s ,+)(d x <j) s ,-) + — ^cos(2<?j s ) 
- 77-77 cos(2(j) c ) cos(2(j) s ). 



(I) 
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Here 4> c / s = 4'c/s.+ + 4 ) c/s,- are the total phase fields, and 
the couplings are given by 

9 P =92± + .92|| - .9i|| = U + 3Ji/2 - 3J 2 /2, (8a) 

9 a =92± - .92|| + .9i|| = U- Ji/2 - 3J 2 /2, (8b) 

,9c = 93± = U + 3J 1 /2- 3 J 2 /2, (8c) 

g s =g lx = U- Jx/2-3J 2 /2, (8d) 

9cs=g 3 \\=-Ji/2 + J2/2. (8e) 

The renormalized Fermi velocities are v p = 2< + (<74|| + 
<?4±)/2t>" and iv = 2t + (g 4 || — 34_l)/27t. Here and in the 
following, we use the lattice constant as the unit for the 
coupling constants as well as for the Fermi velocities. 

The spin-charge coupling term with coupling constant 
g cs describes umklapp scattering processes between elec- 
trons with the same spin. This interaction term formally 
occurs in the zeroth order of the expansion of the fermion 
fields with respect to the lattice constant. It is clear, how- 
ever, that (7 3 || type processes can give contributions only 
for non-local interactions. Moreover, this spin-charge 
coupling term breaks the global spin SU(2) symmetry 
of the system. Therefore, in order to preserve this sym- 
metry, and in order to treat the non-local interactions in 
a consistent way, the next-to-leading terms have to be 
taken into account in the expansion of the fermion fields. 
To first order, among other contributions, three new spin- 
charge coupling terms appear in the Hamiltonian. Wc 
find that the spin and charge velocities are changed by 
the term (— 5 , 1 ||/2), and the symmetry-restoring non-local 
interaction terms are given by 

H (1) (x) =|^(&0 s ,+)(&0 s ,_)cos(20 c ) 

- j^(d x <p Ct+ )(d x <p c ~)coti(2(j) s ) 

+ ^(d x cj> c , + )(d x (l> c ,-){d x cl> s , + ){d x <i> 8 ,-). (9) 

The first two terms correspond to backward and umklapp 
scattering, respectively, between electrons with opposite 
spins, and the third term describes backward-scattering 
processes between electrons with equal spins. Initially, 
all these couplings are equal to g cs , 



9 P s 



= 9pa=9cs = -Ji/2 + J 2 /2. 



(10) 



The SU(2) symmetry of the spin sector assures g s = g ai 
9cs = 9c<r, and g ps = g pa . Therefore, there are five inde- 
pendent couplings which we choose to be g p , g c , g s , g cs , 
and g ps . We note that the renormalization of the Fermi 
velocities, which is a secondary effect, will not be taken 
into account in the following. 



C. Renormalization group analysis for fluctuating 
charge and spin fields 

The Hamiltonian H = + ijW cannot be solved 
exactly However, a renormalization group (RG) analy- 
sis permits the investigation of the relative importance 



of the various couplings. In the RG procedure, the cou- 
plings are considered to be a function of some scaling pa- 
rameter y, e.g., the logarithm of the effective bandwidth. 
As the scaling parameter is taken to infinity, the flow of 
the couplings shows which of them are important and 
which can be ignored, depending on whether or not they 
tend to zero, to a finite value, or to infinity. For exam- 
ple, when all couplings but the forward scattering terms 
tend to zero, the Hamiltonian H describes a Luttingcr 
liquid with freely propagating charge and spin degrees of 
freedom. 

The one-loop RG equations for our five dimensionlcss 
running coupling constants g x (y) = gx{y)/^t rea d 15 i 18 



dffp(^) 

dy 
dffc(y) 

dy 

dg s (y) 

dy 

dg cs (y) 
dy 

&9 P s{y) 

dy 



2gl + & + 9s9 P s, 



■ ^g P gc - 9s9cs - g cs g P 
- 2.9? ~ ~9c9cs - g 2 s , 



(lla) 
(lib) 
(11c) 



2g cs + 2g p g cs - 4g s g cs - 2g c g s 
^9c9 ps - ^9cs9 ps , (Hd) 

^9 PS + 2g P g s - ^gc9cs - ^g 2 cs 

^9s9ps , (He) 



with initial values g x {y = 0) = g x /ATtt. From these equa- 
tions, it follows that there is only a single line of weak- 
coupling fixed points, namely ~g c = Jj s = ~g cs = ~g s = 0. 
In order to show this, we note that we have started our 
analysis assuming that there is neither a charge gap nor 
a spin gap. This implies that a weak-coupling fixed point 
corresponds to g c = g s = 0. Equations (jlip immediately 
imply that ~g cs = ~g = also, and that only ~g p remains 
undetermined. 

A linear stability analysis of the fixed-point line shows 
that it is stable against small perturbations g cs and g ps , 
that it is marginally stable against small perturbations 
g s and g p , and that its stability with respect to pertur- 
bations g c depends on the sign of the fixed-point value 
g p (stable for g p < 0, unstable for g p > 0). Therefore, in 
order to determine the weak-coupling regime, it is conve- 
nient and sufficient to consider the RG equations without 
the spin-charge coupling terms, i.e., we may consider the 
RG equations for g cs = g ps = 0. We thus arrive at 



dffp(y) 
dy 

^■9 civ) 
dy 

dg s (y) 

dy 



igl 



= ^9 P 9c 



= - 2.9.? 



(12a) 
(12b) 
(12c) 



in the vicinity of the weak-coupling fixed-point line. 

This simpler problem is readily analyzed. The trajec- 
tory for the spin coupling g s (y) flows to infinity if g s < 0. 
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In this case, a gap opens in the spin spectrum. If g s > 0, 
this coupling is marginally irrelevant, i.e., the spin mode 
remains soft. In the charge sector, g p = g c initially, and 
this relation remains valid under the RG flow. Therefore, 
it is sufficient to consider Eq. (|12a[) . It is seen that for 
g c > the charge mode becomes gapped because g c (y) 
flows to infinity, otherwise the charge excitations remain 
gapless. 

The simplified equations show that a fully gapless 
Luttinger-liquid phase, ~g c = ~g s = 0, is not possible for 
our model. The initial couplings would have to fulfill 
g c < and g s > which requires J 2 > 2 £7/3 + J\ for 
g c < and J 2 < (2U — Ji)/3 for g s > 0. These two con- 
ditions cannot be fulfilled simultaneously with positive 
bare couplings U, Ji, and J2. Consequently, we must 
redo our RG analysis under the assumption that at least 
one of the two modes is gapped. 



D. Renormalization group analysis for gapped 
charge and/or spin modes 

When one of the fields is gapped, the spin-charge cou- 
pling processes become relevant ! 18 ' 19 Their contribution 
will be considered on the mean-field level. In this picture, 
the gapped field is locked to a value which optimizes the 
interaction energy. 

When there is a gap in the charge sector, the charge 
field 4> c is locked at <fi c = 0mod7r because the initial 
value of the coupling g c is positive. Neglecting the fluc- 
tuations of the field cf> c in Hamiltonian ((9|), the terms pro- 
portional to g pa and g prT do not contribute, and cos(2$> c ) 
can be replaced by its weak-coupling mean-field value, 
cos(20 c ) = 1. Due to this substitution, the interaction 
terms proportional to g cs and g c „ become marginal be- 
cause their scaling dimensions reduce to x cs = x ca = 2. 
On the mean-field level, the spin-coupling term propor- 
tional to g cs is of the same form as the interaction term 
proportional to g s in H (°) . Therefore, the spin field (f> 8 
fluctuates in the modified potential g* cos(20 s ) with the 
new coupling g* s , 



9* 8 = 9s - 9cs = U - 2J 2 . 



(13) 



Analogously, the interaction term proportional to g ca in 
combines with the interaction term proportional to 
g a in H(°> to produce the new coupling g*, with 



9a = 9a 



U-2J 2 



(14) 



This equation shows that the SU(2) spin symmetry is 
preserved on the mean-field level. 

In the presence of a charge gap and the SU(2) spin 
symmetry, we only have to analyze a single equation for 
g s instead of the five RG equations (fTTj) . namely 



dy 



-25. 



S 1 



(15) 



with the initial value g s (y = 0) = g*/4iTt. It is readily 
seen that the spin mode becomes gapped if g* < 0, i.e., 



J2 > U/2, independently of the value of the nearest- 
neighbor interaction J\. 

When there is a gap in the spin sector, the spin field 
(j) s is locked at <j) s = mod tt because the initial value of 
the coupling g s is negative. Neglecting the fluctuations 
of the field <fi s in the Hamiltonian ©, the terms pro- 
portional to g pa and g ca do not contribute and cos(20 s ) 
can be substituted by its weak-coupling mean-field value, 
cos(20 s ) = 1. Due to this substitution, the interaction 
terms proportional to g cs and g ps become marginal be- 
cause their scaling dimensions reduce to x cs = x ps = 2. 
On the mean-field level, the charge-coupling term pro- 
portional to g cs is of the same form as the interaction 
term proportional to g c in . Therefore, the charge 
field 4> c fluctuates in the modified potential g* cos(2(/> c ) 
with the new coupling g* c , 



9 c 



9c + 9cs = U + Ji - J 2 



(16) 



Using similar reasoning, the new coupling g* becomes 

g* p = g P -g ps = U + 2J 1 -2J 2 . (17) 

Note that these new initial couplings are not equal, so we 
must analyze the two-dimensional scaling curves defined 
by the equations 



dy 

<jg c {y) 

dy 



2gl 



2g P g c , 



(18a) 
(18b) 



given the initial values g c (y = 0) = g*/4irt and g p (y = 
0) = <7*/47rf. The flow diagram is shown in Fig. [TJ 




FIG. 1: Scaling curves for the charge-coupling parameters g c 
and g p in the presence of a spin gap. 

The conditions for a gapped charge mode are either 
g* > or (g* < and \g*\ > \g*\). This leads to the 
result that a gapped charge mode exists if J2 < 2U/3+ J\. 



E. Phase diagram 

In general, we find three regions where either the 
charge gap or the spin gap or both are finite. It is inter- 
esting to analyze the dominant correlations in the var- 
ious gapped phases. The order parameters for density 
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waves in the charge (CDW), spin (SDW), bond-charge 
(BCDW), and bond-spin (BSDW) require the calcula- 
tion of correlation functions using the operators 

0i,CDW = (-I)'Kt + m,l), (19a) 
Oi,SDW = (-l) i K T -n u ), (19b) 

C\bcdw = (-l)*(ct + c] d c i+lt i + h.c), (19c) 

Oj.BSDW = (-l)*(c],1-Ci+l,T - cT^Cj+1,4, + /i.e.) , (19d) 

written in terms of the lattice fermions. These order 
parameters become 

Ccdw(^) oc sin 0c (x) cos</> s (x), (20a) 

Osdw(x) oc cos</> c (x) sin</> s (x), (20b) 

Cbcdw {%) oc cos <f> c (x) cos 4> s (x) , (20c) 

Obsdw (x) oc sin C (x) sin </> s (x) (20d) 

in bosonized form. When the charge mode is gapped, the 
field (f> c is locked at cf> c = mod ir. When the spin mode is 
gapped, the field <p s is locked at <fi s = mod ir. Therefore, 
in the regime where both of the fields are gapped, we find 
that the BCDW order parameter is maximal. Therefore, 
the model describes a phase with bond ordering (BOW) 
for A c ^ and A s ^ 0. 

When only the charge mode is gapped, the spin field is 
a free field. However, upon increasing the scaling param- 
eter (y) of the renormalization group procedure, the ini- 
tially negative spin coupling grows and tends to zero, and 
the spin field oscillates around ir/2 (mod ir). Therefore, 
for small couplings, the dominating ordering is SDW for 
A c 7^ and A s = 0. Note that the SU(2) spin symmetry 
is not spontaneously broken, i.e., the spin correlations 
are critical without true long-range order. 

Similarly, when the spin mode is gapped and the charge 
mode is gapless, there is no true long-range charge or- 
der. Therefore, we call this phase the Luther-Emery (LE) 
phase. The charge coupling g c tends to zero, either from 
positive values or from negative values. Depending on 
the sign of the charge coupling, <p c fluctuates around 7r/2 
or around zero. Correspondingly, the dominating correla- 
tions are either CDW or BCDW for A c = and A s ^ 0. 
The line which separates the dominant BCDW critical 
correlation and the dominant CDW correlations in the 
LE phase is indicated in Fig. [5] by a dashed line. 

The resulting phase diagram of the t-U-J\-Ji model at 
weak coupling is shown in Fig. [2j For U = 0, the spin 
gap is always finite for Ji > 0. For Ji < Ji, the charge 
gap is also finite, and the ground state is characterized 
by a bond-order wave. The charge gap closes at Ji = J\ 
and the system goes into a LE phase with no long-range 
charge or spin ordering but critical charge-density-wave 
correlations. 

For U > 0, Ji > 0, and Ji < U/2, the ground state 
is analogous to the spin-density-wave (SDW) phase of 
the one-dimensional Hubbard model, i.e., the charge gap 
is finite, the spin gap is zero, and the spin correlations 
are critical. For 2U/3 + J± > J2 > U/2, both the spin 
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A c ?0 
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FIG. 2: Field-theory prediction for the half-filled t-U-Jx- 
J'2 model. The solid lines give the phase boundaries be- 
tween the fully gapped regime (bond-order wave, BOW) and 
the semi-gapped regimes (spin-density wave, SDW; Luther- 
Emery, LE) . The dashed line shows the border between dom- 
inantly charge-density-wave and bond-order-wave correlations 
in the Luther-Emery phase. 



gap and the charge gap are finite. The ground state is a 
BOW with long-range order in the bond-charge-density- 
wave correlations. For J2 > 2L//3 + Ji, the charge gap 
closes and the system goes over to the LE phase with 
a finite spin gap but no charge long-range order. For 
2J7/3 + J\ < J2 < U + Ji, the bond-charge-density-wave 
fluctuations dominate, whereas, for J2 > U + Ji, the 
fluctuations in the charge-density-wave order parameter 
are maximal. 



u 



SDW 



A c ^0 
A s = 



FIG. 3: Field-theory prediction for the half-filled t-U-Ji 
model. For all Ji > 0, the ground state is a spin-density- 
wave (SDW) phase with a finite charge gap, zero spin gap 
and critical spin correlations. 



In order to make contact with earlier work, we dis- 
play the phase diagram of the t-U-J\ model separately 
in Fig. [3J In contrast to previous resultS ) 14 ' 15 i 16 we do 
not find any signature of a BOW phase. For all J\ > 0, 
the ground state is SDW, just as is the ground state of 
the half-filled Hubbard model for U > 0. This result is 
corroborated by our numerical DMRG data, which we 
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present in the next section. 



III. NUMERICAL RESULTS 

In order to explore the phase diagram of the Hamilto- 
nian ([1]) and to test the predictions of bosonization, we 
carry out extensive, high-precision, ground-state DMRG 
calculations. 20 ' 21 ' 22 Relatively high sensitivity is required 
to resolve the phases, especially in the weak-coupling 
regimes in which one would expect bosonization to be 
valid. In order to differentiate the possible phases, we 
calculate the spin gap A s , the charge gap A c , and the 
bond-order-wave parameter (B) of the one-dimensional 
t-U-J\-J2 model on lattices with open boundary condi- 
tions and up to L = 256 sites. The weight of the dis- 
carded density-matrix eigenstates is held below a maxi- 
mum of 10 -9 . 

For finite systems, the spin gap A S (L) is defined as 

A S {L) = E (L, N,S = 1) — Eo(L, N,S = 0). (21) 

Accordingly, the charge gap A C (L) is determined using 

A C (L) = [E Q (L,N + 2,S=1)+E (L,N-2,S = 0) 
-2E (L,N,S = 0)}/2, (22) 

where Eq(L, N, S) is the ground-state energy for an L- 
site system with N electrons and total spin S. We ex- 
trapolate using second-order polynomials in l/L to de- 
termine the spin gap A s and the charge gap A c in the 
thermodynamic limit, 
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FIG. 4: (Color online) Finite-size extrapolation of the spin 
gap as a function of l/L for the t-U-Ji model at (a) U = 
and (b) U = 0.1. 



Results for J2 = 



A S (L) = Af +A S /L + B S /L 2 1 

A C (L) = A™ +AJL + BJL 2 , (23) 

where A£° s , A CjS , and B c s are fitting parameters. The 
staggered bond order parameter is defined as 

(B)(L) = _L^^^(-1)'+ 1 (4 Cj+1iCT+ /,. c ). (24) 

The bond order parameter (B) is extrapolated using 
finite-size corrections of the form l/L 1 1 without consid- 
ering higher corrections, 

(B)(L) = (B°°)+A B /L\ (25) 

where (B°°), As, and 7 are fitting parameters. We find 
that adding higher-order terms, which increases the num- 
ber of fit parameters, tends to make the fits less stable. 

In the following, we first treat the t-U-Ji model, i.e., 
J2 = in Hamiltonian ([1]), then study finite positive J2, 
first with [7 = 0, then with nonzero U . For simplicity, 
in the remainder of this article the energy scale is set by 
taking t = 1, and so U, J\, and J2 are dimcnsionlcss 
quantities. 



For the unfrustrated case (J2 = 0), our bosonization 
procedure of Sec. [TT] predicts a SDW phase with a finite 
charge gap and critical gapless spin excitations, A c > 
and A s = 0. In the SDW phase, the bond order param- 
eter vanishes. 

The finite-size extrapolation of the spin gap, plotted 
as a function of l/L for U = and U = 0.1 is shown 
in Fig. [4] As can be clearly seen, the scaling behavior is 
predominantly linear in l/L, and the l/L — > extrapo- 
lated value, A^°, is zero on the scale of the plot for all 
values of Ji for both values of U. A fit of the data with 
a second-order polynomial in l/L, as discussed above, 
yields a value of A s that is less than 2 x 10 -4 in all cases. 
This puts a rather stringent constraint on bond ordering 
in this case; the spin excitations are gapless to a very 
high numerical accuracy. 

The system-size behavior of the charge gap is displayed 
in Fig. [5l As can be seen, the l/L — > extrapolated 
value, A£°, is nonzero in general, with the scaling going 
from being predominantly linear in l/L (with a small 
negative (l/£) 2 term) when A c is small, to having a 
substantial positive (l/L) 2 term when A c is significantly 
different from zero. Such finite-size behavior is typical 
for gaps in one-dimensional systems with open boundary 
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FIG. 5: (Color online) Finite-size extrapolation of the charge 
gap as a function of 1/L for the t-U-Ji model at U = 0. 
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FIG. 6: (Color online) Extrapolated spin and charge gaps for 
the t-U-Ji model at U = and 0.1 as functions of Ji. 



ciablc positive finite value which varies significantly as a 
function of J\. The L — > oo extrapolated value is 
small, but still shows some variation with Ji . Note, how- 
ever, that the extrapolated value is negative at small and 
large J\ and is positive only for intermediate J\. Taking 
the largest negative value ((-B°°) ~ —0.003) as a rough 
estimate of the extrapolation error, the largest positive 
value, (B°°) « 0.007, is not distinguishable from zero to 
within our accuracy. Moreover the fit to Eq. (I25[) yields 
an exponent 7 which varies between 0.47 and 0.77. All 
this underlines the uncertainty in carrying out extrap- 
olations using this analytic form and the sensitivity of 
(B°°) to the details of the fit. On the other hand, as 
discussed above, vanishes to a high accuracy for all 
Ji, precluding a BOW phase. Thus, within the numer- 
ical methods applied here, the spin gap seems to be a 
significantly more sensitive probe for the existence of a 
bond order wave phase than the bond order parameter 
(B) itself. 
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conditions. 

The behavior of the extrapolated gaps as a function of 
Ji is shown in Fig. [6] As discussed above, the spin gap 
is numerically indistinguishable from zero for all values 
of Ji for both U — and U = 0.1. The extrapolated 
charge gap is small on the scale of the plot for J% < 0.8, 
and then increases, crossing over to a linear increase for 
larger values of J\ . From bosonization, we would expect 
an exponential opening of the gap with J\, similar to 
the exponential opening of the charge gap with U in the 
Ji = c&scM The Ji-dependence of in Fig. [6] is 
qualitatively consistent with such a behavior. We have 
not carried out an explicit fit because the detailed form of 
the exponential opening is not known from bosonization; 
to determine the specifics of a general exponential form 
via fitting to finite-size extrapolated data is difficult. 

We now turn to the BOW order parameter, displayed 
as a function of J\ for various system sizes and L = 00 
in Fig. [7J At each system size, (B)(L) has an appre- 



FIG. 7: (Color online) Bond-order parameter (B)(L) for 
L — 32, 64, 96, 128, 256, 512 and extrapolated bond order pa- 
rameter {B°°) as a function of Ji for the t-U-Ji model at 
(7 = 0. 



Our DMRG calculations for J2 = are thus in agree- 
ment with the predictions of the bosonization calcula- 
tions of Sec. [Til sec Fig. [3J the ground-state phase is a 
SDW with gapless spin excitations for all positive U and 
J\. While we have treated explicitly only two values of 
the interaction strength, U = and U = 0.1, we have 
chosen these values in accordance with the phase dia- 
grams of Refs. [H and Qjl which predict the appearance 
of a bond order wave phase only for U < 0.35. At larger 
values of U, the behavior should be that of the ordinary 
half-filled Hubbard chain and one would not expect a 
BOW phase to occur. 
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FIG. 8: (Color online) Finite-size scaling analysis for (B)(L) 
for different Ji when U — and Ji = 1. The inset shows 
the finite-size scaling analysis for J2 = 2.0, 2.5, and 3.0, from 
bottom to top. 
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FIG. 10: (Color online) The L = 00 extrapolated bond order 
parameter {B°°) as a function of J2 for U — and Ji = 1. 
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FIG. 9: (Color online) Extrapolated spin gap and charge gap 
as functions of J2 for U = 0, Ji = 1. The inset displays the 
same data for J2 < 2.5 on an enlarged scale. 



B. Results for (7 = and nonzero J2 

We now include the explicit frustration J 2 while setting 
the on-site Coulomb interaction to zero. Fig.[5]shows the 
system-size extrapolated spin and charge gaps, and 
AJ?°, as functions of J 2 at U = and J\ = 1. (We 
do not show the finite-size extrapolation, which proceeds 
similarly to that in Figs. 2] and [H explicitly.) The spin 
gap opens slowly at small J 2 , but with a form consistent 
with a critical Jf = (see the inset in particular). The 
charge gap decreases rapidly with J 2 at small J 2 , reaching 
zero at J% ks 1 = J\, but then opens again at J 2 ^ ~ 2. 
At weak coupling, this behavior of both gaps is consistent 
with the predictions of bosonization, but the reopening 
of the charge gap for larger J 2 is not contained in the 
bosonization analysis. However, such large values of J 2 
are clearly outside its region of validity. 



Representative results for the finite-size scaling of the 
bond order parameter (B) are present in Fig. [HJ For 
small J2, the scaling behavior is similar to that for J 2 = 
0, yielding an exponent 7 that varies between 0.44 and 
0.71. However, for large J 2 , the data extrapolate almost 
linearly to finite values. This illustrates that the scaling 
form (|25|) . goes over to a function that might be better 
fit by a polynomial in 1/L, as in Eq. (p?5|) . However, for 
consistency, we nevertheless always use Eq. (f2~5j) for the 
fitting and note that the case of a linear function of 1/L 
is encompassed by Eq. with 7 = 1. 

The extrapolated results for (B°°), plotted as a func- 
tion of J 2 , are shown in Fig. |TD] For J 2 = to J 2 ^\ 
(B°°) is very small, even falling off from the small finite 
value at J 2 = 0, which we have argued to come about due 
to numerical and extrapolation errors. Note that here, 

c(l) 

for J 2 < J 2 ' ?s 1, the phase is characterized as bond 
order wave within bosonization. While this seems to be 
a contradiction at first glance, note that the charge gap, 
Fig. O falls off very rapidly from its small finite value at 
J 2 = 0, whereas the spin gap opens very slowly due to 
its putative exponential form. In consequence, the value 
of (B°°) is very small. Our interpretation, then, is that 
the BOW order parameter is finite, but numerically un- 

c(l) c(2) 

resolvable in this region. For J 2 y± > <J 2 < J 2 V ', the spin 
gap is clearly non- vanishing, but (B°°) is numerically 
zero. This behavior is consistent with the bosonization 
prediction of a Luther-Emery phase. In other words, the 
vanishing charge gap indicates a phase in which there is 

c(2) 

no BOW. When J 2 > J 2 , coincident with the reopen- 
ing of the charge gap in Fig. [51 the BOW phase reappears, 
this time clearly marked by a finite bond order parameter 
as well as finite spin and charge gaps. 
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FIG. 11: (Color online) The L = oo extrapolated spin gap FIG - 12: ( Color online ) The L = oo extrapolated (B) as a 
and charge gap as functions of J 2 for U = 2, Ji = 1. function of J 2 for U = 2, Ji = 1. 



C. Results for nonzero U and J2 

We now study the effect of the frustration J 2 when the 
Coulomb repulsion U is finite. Bosonization predicts that 
the SDW phase that is present only along the J2 = line 
at U = becomes enlarged to a finite region at finite U. 
We explore the behavior as a function of J 2 for moderate 
values of U and J 1: U = 2, and Ji = 1. Fig. [TT1 shows 
the spin and charge gaps, extrapolated to infinite systems 
size, as a function of J 2 . As can be seen, the spin gap 
opens at a finite J 2 « 0.6 and the charge gap, although 
at first decreasing and reaching a minimum at J2 ~ 1.1, 
is always finite. As can be seen in Fig.[T2J the bond order 
parameter (B°°) = when J2 < Jf , and opens rapidly to 
a large, finite value at J 2 ~ 0.5. The behavior of all quan- 
tities is consistent with a SDW phase for small J2 and 
a BOW for large J 2 . Bosonization does predict a tran- 
sition from a SDW phase to a BOW phase at J2 = U /2 
(see Fig. [2]). However, it also predicts a transition to a 
spin-gapless LE phase at larger J2, which is not found 
in the numerical calculations. In our opinion, this is be- 
cause the values of U, J\ and J 2 here arc large enough so 
that the regime of validity of bosonization is exceeded. 
Note that the critical value JJ sw 0.5 is far from the 
weak-coupling prediction of J2 = U/2 = 2, but agrees 
fairly well with the value expected from the frustrated 
Heisenberg chain, for which {Jf eis / Jf eis ) c « 0.241^2^ 
if we take Jf eis = J x + At 2 /U = 3, the effective Heisen- 
berg coupling within strong coupling; this yields an esti- 



mate J. 



c(strong) 



0.72, in reasonable agreement with the 



DMRG result. 

Fig. summarizes the phase diagrams as a function 
of J2 obtained from the DMRG calculations at zero and 
finite U. For {7 = 0, the SDW phase at J 2 = be- 
comes a BOW phase at arbitrarily small, but weak J%. 
At intermediate J2, a metallic, but spin-gapped Luther- 
Emery phase occurs, and at large J2 the system reenters 
the BOW phase. At moderate, finite U, the SDW phase 
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FIG. 13: A sketch of the ground-state phase diagram of the 
one- dimensional t-U-.J1-.J2 model at zero and finite U ob- 
tained from analysis of the DMRG calculations. 



persists when J2 is small and finite, going over to a BOW 
at larger J 2 . 



IV. DISCUSSION AND CONCLUSION 

In this work, we have investigated the ground-state 
behavior of the half-filled one-dimensional Hubbard 
model with antifcrromagnctic nearest-neighbor and next- 
nearest-neighbor Heisenberg interactions. Our field- 
theoretical analysis for weak couplings indicates that the 
ground state has a finite gap for either charge excita- 
tions (spin-density-wave phase, SDW) or spin excita- 
tions (Luther-Emery phase, LE) or both (bond-order- 
wave phase, BOW). Our extensive numerical DMRG in- 
vestigations agree very well with the field-theoretical pre- 
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dictions for small interactions. The only exception is the 
lack of numerical evidence for a finite bond-order param- 
eter in the region U = 0, Ji = 1 and < J2 < Ji. Here 
the system sizes are large enough to resolve finite spin 
and charge gaps but they are still too small to detect the 
very small bond order parameter. 

For larger interactions, e.g., U = 2, the DMRG finds a 
strong-coupling bond-ordcr-wave phase which eludes the 
field-theoretical description. Instead, its existence and 
its properties can be inferred from a strong-coupling ex- 
pansion of the model where it is seen that the strong- 
coupling BOW phase results from the frustration of the 
nearest-neighbor and next-nearest-neighbor Heisenberg 
couplings. Therefore, the metallic Luther-Emery phase 
is limited to a narrow weak-coupling region in the phase 
space where it would be very difficult to justify the 
strengths of the coupling parameters from microscopic 
considerations. For moderate interactions, an echo of 
the weak-coupling Luther-Emery phase can be seen in 
the behavior of the charge gap as a function of J2 , which 
displays a minimum at some Ji > J\. 

The nearest-neighbor Heisenberg coupling J\ is not a 
frustrating interaction for the half-filled Hubbard model 
because the ground state of the t-U-J\ model is a spin- 
density wave for all Ji > 0. In order to arrive at this 
conclusion in the field-theoretical analysis, the fact that 
bosonic phase fields are locked to their mean-field val- 
ues when excitations are gapped, so that seemingly ir- 



relevant operators become marginal operators, must be 
taken into account. In numerical calculations one needs 
to study rather large system sizes in order to extrapo- 
late to a vanishing spin gap and bond-order parameter 
in the thermodynamic limit. The next-nearest-neighbor 
Heisenberg interaction Ji, in contrast, truly frustrates 
the Hubbard model, opening the way to Luther-Emery 
and bond-charge-ordcred phases for J2 > 0. As ex- 
pected from our experience with the frustrated Heisen- 
berg model, the SDW phase is stable against weak frus- 
tration for U > 0, i.e., a finite J2 is required to open the 
spin gap. 

In conclusion, our study demonstrates both analyt- 
ically and numerically that a nearest-neighbor Heisen- 
berg exchange interaction added to the half-filled Hub- 
bard model does not lead to frustration or to new phases 
in the ground-state phase diagram, whereas a frustrating 
next-nearest-neighbor exchange docs. 
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